The X X S + and a 3 S + states of LiCs studied by Fourier-transform spectroscopy 
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We present the first high-resolution spectroscopic study of LiCs. LiCs is formed in a heat pipe oven 
and studied via laser-induced fluorescence Fourier-transform spectroscopy. By exciting molecules 
through the X 1 E + -B 1 n and X 1 S + -D 1 II transitions vibrational levels of the X 1 E + ground state have 
been observed up to 3 cm -1 below the dissociation limit enabling an accurate construction of the 
potential. Furthermore, rovibrational levels in the a 3 E + triplet ground state have been observed 
because the excited states obtain sufficient triplet character at the corresponding excited atomic 
asymptote. With the help of coupled channels calculations accurate singlet and triplet ground state 
potentials were derived reaching the atomic ground state asymptote and allowing first predictions 
of cold collision properties of Li + Cs pairs. 

PACS numbers: 31.50.Bc, 33.20.Kf, 33.20.Vq, 33.50.Dq 



I. INTRODUCTION 

Spectroscopy of heteronuclear diatomic alkali 
molecules provides important input to current re- 
search in cold molecules and mixtures of ultracold 
atomic gases. Cold heteronuclear alkali dimers are 
subject to a large interest since they can be formed at 
temperatures below 1 mK 0, [H, |j, || and possess a 
large permanent electric dipole moment for deeply bound 
singlet levels @, 0| • This combination of properties en- 
ables electric field control of ultracold collisions [1, 
and cold chemical reactions fiol fill ] and holds promises 
for applications in quantum computation fl2l |. Precise 
potential energy curves, in particular for the lowest 
electronic states, are evidently important for such 
applications as well as for understanding the molecule 
formation processes (photoassociation), ro- vibrational 
state selective detection [f| [l3[ and for formation of 
vibrational ground-state molecules [3]. In ultracold 
mixtures of atomic gases, quantum degeneracy of one 
atomic species can be achieved through sympathetic 
cooling by the other species [HI, [13] • The interspecies 
interaction strength can be varied through magnetic 
Feshbach resonances and has a large influence on a 
variety of effects such as phase-separation between a 
Bose-Einstein condensate and a degenerate Fermi gas 
[TH and the transition to a Bardeen-Cooper-Schrieffer 
superfluid state in dilute Fermi gases Understand- 
ing interspecies collision properties at the atomic ground 
state asymptotes, e.g. predicting the magnetic field 
strength values of Feshbach resonances at zero kinetic 
energy, requires precise potential curves of the electronic 
ground states. 

Cold LiCs molecules were observed very recently, being 
formed in a two-species magneto-optical trap Previ- 
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ously, inelastic collisions [20| and sympathetic cooling of 
Li by Cs in an optical dipole trap [2l| have been studied. 
Recent theoretical work considers LiCs in strong dc elec- 
tric fields and its influence on rovibrational dynamics [22| 
and Li-Cs collision cross sections [9j] . 

Although LiCs was observed already in 1928 by Walter 
and Barratt through absorption in a mixture of metal- 
lic vapors (23| . very few and no high-resolution spectro- 
scopic studies have been made until now. In the 1980's 
Vadla et al. studied the repulsive Li(2p)+Cs(6s) asymp- 
tote [24j]. More recently LiCs molecules were formed on 
He nanodroplets and excitation spectra of the d 3 n <— 
a 3 S + transition were recorded and modeled (25[. Ab ini- 
tio potentials were calculated by Korek et al. [26[ (see 
Fig. [T]) and more recently an extended theoretical study 
was done by Aymar and Dulieu The theoretical po- 
tentials provide a good starting point for analysis of the 
spectra obtained in the present work. 

Here we present a high-resolution spectroscopic study 
of the LiCs molecule. Similar to our previous studies 
[UHl] we apply Fourier-transform spectroscopy of laser- 
induced fluorescence from LiCs molecules formed in a 
heat-pipe, because this technique is suitable for collect- 
ing a large amount of accurate experimental data. By us- 
ing properly chosen excitation schemes (see e.g. (27l. [28j) 
we can measure transition frequencies to a wide range 
of vibrational levels in the singlet X 1 !]" 1 " as well as in 
the triplet a 3 S + ground states, especially levels close to 
the asymptote. Since the two states are coupled at long 
internuclear distances by the hyperfine interaction it is 
not correct to treat them separately in this region, as it 
would lead to model potentials which are unable to repro- 
duce the experimental observations close to the asymp- 
tote. Therefore, the aim of our experimental work is 
to collect experimental data on both ground states and 
to fit accurate experimental potential energy curves si- 
multaneously for both states - the indispensable starting 
point for a study of the molecular structure of LiCs or for 
modeling of cold collision processes on the Li(2s)+Cs(6s) 
asymptote. We apply these potential curves within a cou- 
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FIG. 1: (Color online) Selected singlet and triplet potential 
energy curves in LiCs [26fl . 

pled channels model in order to explain our experimental 
observations and also to compute collision properties for 
comparing recent results of sympathetic cooling of Li by 
Cs [lU 

The article is organized as follows: The experimental 
setup and excitation schemes are presented in Section HT1 
In Section Mil we describe the analysis of the obtained 
spectra. The procedure for construction of potential en- 
ergy curves is described in Section HVl and the potentials 
are reported. In Section [V] we give our conclusion and 
an outlook for further experimental study needed for a 
quantitative description of ultracold collisions in Li + Cs. 



II. EXPERIMENT: MOLECULE FORMATION 
AND SPECTROSCOPY 

A. Molecule formation 

LiCs molecules are formed in a stainless steel heat-pipe 
identical to the one described in Ref . [29} , except for one 
modification described below. The heat-pipe (960 mm 
long and 34 mm outer diameter) is filled with 6 g Li and 
5g Cs (the CSjhi a closed ampoule, is loaded into the 
side container [29| Q) and typically operated with 3-6 
mbar Ar buffer gas pressure. 

Since the vapor pressure difference between Li and Cs 



at a common temperature is extremely large, we mod- 
ified the design of Ref. 29] in order to obtain a three- 
section heat-pipe, which is more suitable for producing a 
vapor mixture with similar concentrations of Li and Cs 
and hence for forming LiCs molecules [3(j. In Ref. [29[ 
the central 60 cm of the heat-pipe is heated uniformly 
in a commercial oven (Carbolite); here we mount two 
stainless steel 'shells' (20 cm long, 50 mm inner diameter) 
concentrically around the heat-pipe and seal up the ends 
facing towards the center of the oven such that only the 
central 20 cm of the heat-pipe are heated directly in the 
oven. By blowing air into the open ends, which extend 
outside of the oven, we can maintain a lower temperature 
in the sections shielded by the shells than in the central 
part of the heat-pipe. 

The heat-pipe is conditioned by heating it to temper- 
atures of about 580 °C under 10 mbar Ar pressure; sub- 
sequently the Cs ampoule is broken by shaking the tube. 
Operating temperatures are 540 °C in the central part 
and 370 °C in the outer sections. The heat-pipe oven 
was operated for more than 200 hours over a 10 month 
period and was still in good working conditions at the 
end of this period. 



B. Laser-induced fluorescence Fourier-transform 
spectroscopy 

Laser-induced fluorescence from LiCs molecules is ob- 
served after excitation on the B 1 !! <— X 1 S + and D 1 !! <— 
X X X+ transitions. The B 1 !! <- X 1 S+ transitions were 
excited using a Coherent 599 dye laser (with DCM dye) 
at frequencies in the range 15529-16123 cm -1 and a Co- 
herent 699 dye laser at frequencies in the range 16397- 
17022 cm" 1 (with Rhodamine 6G dye as well as with 
a mixture of Rhodamine 6G and Rhodamine B). The 
D 1 !! <— X : I] + transition was excited using the dye laser 
with Rhodamine 6G at frequencies in the range 16663- 
17238 cm -1 and a frequency doubled Nd:YAG laser. The 
strongest signals were observed for the B 1 !! <— X 1 S+ sys- 
tem and we studied the ground states mainly through 
this system. We note that indeed Walter and Bar- 
ratt observed strong absorption in the range 15983- 
16582 cm -1 [23| which according to Fig. Q] corresponds 
to the B 1 !! <— X 1 E + transition. At excitation frequen- 
cies in the range ~ 14900 — 15400 cm -1 , we find strong 
fluorescence due to Li2 and NaCs (Na is present as an 
impurity in the Li sample) which overshadows a possi- 
ble LiCs signal. Using an Ar-ion laser for excitation at 
457.9 nm, 476.5 nm, 488.0 nm, 496.5 nm and 514.5 nm we 
did not observe any fluorescence from LiCs molecules, 
only from Li2, LiNa and NaCs. 

Contrary to the previously studied molecules 
NaRb [Ig and NaCs H, the low lying levels of 
the excited B 1 !! state in LiCs turned out to be almost 
free of local perturbations by the neighboring triplet 
states b 3 II and c 3 E + , which could be expected from 
the theoretical potential curves (see Fig. [T]). As a 
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consequence we were unable to register any transition 
to the triplet ground state from the low lying B 1 ]! state 
levels. We searched instead for access to the triplet 
manifold through local perturbations in the D 1 !! and 
C 1 E + states. Here we used a Coumarine 6 dye laser 
with a typical power of 25 mW. Unfortunately within 
the searched excitation frequency region of 18446 - 19039 
cm -1 we were also not able to register transitions to the 
triplet ground state. 

The only fluorescence to the a 3 £ + state is observed af- 
ter excitation to high lying levels in the B 1 ]! state. These 
transitions are attributed to the long-range change over 
of coupling case for the B 1 !! state itself rather than to 
local mixing of the B 1 !! state at long range with the 
neighboring triplet states [27|, |28| . This conclusion is sup- 
ported first, by the observation that high lying B state 
levels seem to be locally unperturbed. Second, the in- 
tensity distribution of progressions from the high lying B 
state levels to the ground singlet state could be explained 
satisfactory by the Franck-Condon factors between the 
X and the B states, inclu ding the highest v'x levels (con- 
trary to the case of NaRb j27|). So transitions from these 
B state levels to high lying a 3 £ + state levels will become 
also probable when the B state changes its character from 
Hund's case (a) 1 II to Hund's case (c) f2 = 1 at long in- 
ternuclear distances. Indeed, in our spectra we find tran- 
sitions mainly to high u" and none to the bottom of the 
triplet ground state. 

The laser-induced fluorescence light is collected in the 
direction opposite to the one of laser beam propagation 
and recorded by a Bruker IFS 120HR Fourier- Transform 
Spectrometer (FTS). For detection we use a photomul- 
tiplier (Hamamatsu R928) or a Si-photodiode. In order 
to avoid illumination of the detector by the He-Ne laser 
(632.8 nm), used in the FTS for calibration and stepping 
control, a notch filter (8 nm full width at half maximum) 
is introduced in the beam path, which suppresses also the 
fluorescence induced in the corresponding spectral region. 
The resolution of the FTS is typically set to 0.03 - 0.05 
cm . The uncertainty of the line positions is estimated 
to be 1/10 of the resolution. For lines with signal-to-noise 
ratio less than 3 the uncertainty is gradually increased. 
Each spectrum results typically from an average of 10-20 
scans, but the number of scans is varied from 5 to 350 
depending on the signal strength for the features of in- 
terest in the spectra. For improving the signal-to-noise 
ratio, the spectral window for some spectra is limited by 
using color glass filters or interference filters. In order to 
facilitate the identification in such cases we recorded also 
a spectrum at the same excitation frequency without fil- 
ters, but with lower resolution (0.1 cm -1 ) and a smaller 
number of scans. 

A list of all excitation frequencies used in these ex- 
periments together with the assignment of the excitation 
transitions are given in Tables 3 and 4 of the supplemen- 
tary materials [3l|. 



III. ANALYSIS OF SPECTRA 

A. The X 1 E+ state 

Assignment of the large number of transitions, about 
6600, to the X 1 E + state is done in an iterative process; 
by gradually improving the potential for the X 1 !]" 1 " state, 
more and more transitions can be correctly assigned. 

Initially, we identify several strong doublet series which 
are clearly recognized as P-R components from the reg- 
ular development of the doublet spacing. Among these 
series we choose those with similar spacing, i.e., with sim- 
ilar rotational quantum numbers. Based on the theoret- 
ical potential for the X 1 S + state [26[ we make an ini- 
tial guess for the vibrational and the rotational quantum 
numbers of these fluorescence progressions. The quan- 
tum numbers that give the closest agreement with the 
theoretical vibrational and rotational spacings are as- 
signed and a small set of Dunham coefficients is fitted. 
If the fit is successful we can use the fitted coefficients 
to assign new progressions; if not, a reassignment of the 
quantum numbers must be made. After some iterations 
we obtain a self consistent set of assigned experimental 
progressions which can be satisfactory described by a few 
Dunham coefficients. This is a first hint of a correct ro- 
tational numbering. Here we point out the very good 
quality of the theoretical calculations @, [26J; the theo- 
retical rotational numbering needed a correction by only 
two or three units. 

When the list of assigned transitions reaches several 
hundreds we perform the first potential fits. Initially, we 
start with a pointwise potential (defined in Section IIV|) 
based on the theoretical curve 12611 and improve it using 
the procedure described in Ref. [32[ . In the further anal- 
ysis of the spectra we apply this pointwise potential curve 
since outside the range of the fitted v" and J" , it usually 
possesses better predictive properties than the Dunham 
type coefficients. Moreover our experience from previous 
studies [27|, l29( shows that a potential curve which fits 
long vibrational progressions recorded at high precision 
in a wide range of rotational quantum numbers indicates 
the correctness also of the vibrational numbering. 

Finally, the established vibrational numbering was 
confirmed by assigning several progressions for the less 
abundant 6 Li 133 Cs molecule. These progressions fit to 
the experimental potential based only on 7 Li 133 Cs data 
when the appropriate reduced mass is applied [33[ . 

In order to describe the important long-range part of 
the potential, it is of great value to collect data with tran- 
sitions to high-lying levels of the ground state. Since we 
noticed during the measurements that transitions to such 
high-lying levels originate from high-lying vibrational lev- 
els of the B 1 n state, we studied this state in order to 
optimize the experimental conditions for observation of 
near asymptotic levels in the X 1 !]" 1 " state (a report on 
the B 1 ]! state is in preparation 34]). From the avail- 
able experimental term energies of the B x n state levels a 
preliminary potential for the B 1 !! state was fitted. This 
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pled channels calculation for J" — N" = 17 in which 
the hyperfine interaction between the states is included. 
The progression to the X 1 E + state is formed by Q-lincs 
whereas the progression to the a 3 E + state consists of 
transitions to A" =15, 17 and 19. Lines not marked 
in the figure belong to another progression also assigned 
and used in our analysis. 

The rotational assignment of the triplet lines is 
straightforward since we always find the progression to 
the X 1 !]" 1 " state which shares the excited level in the B 1 ]! 
state with the progression to the a 3 E + state. The vibra- 
tional assignment of the transitions to the triplet ground 
state is done in the same way as for the X 1 E + state, 
however, no transitions in 6 Li 133 Cs were observed and 
hence the vibrational assignment relies only on the inter- 
nal consistency of the total procedure. 



FIG. 2: (Color online) Hyperfine structure of a transition to 
the v'a = 9, N'a = 17 level of the a 3 E+ state. The vertical full 
lines indicate the prediction of the coupled channels calcula- 
tion. The dashed vertical line indicates the position of the 
hyperfine-structure free level, which is shifted by 0.035 cm" 
from the central component of the structure. 



potential was then used to predict transition frequencies 
for excitation transitions with large Franck-Condon fac- 
tors. In this way we recorded systematically transitions 
from vf = 24 and 25 to the X 1 E + state for a wide range of 
J', thus adding several ro- vibrational levels with — 49 
and 50 to the ground state dataset. In most cases such 
high- lying levels of the X : E + state could be observed 
with a sufficient signal-to-noise ratio only by increasing 
the number of summed scans to several hundreds. We 
searched for excitations to higher vibrational levels of 
the B 1 !! state within predicted spectral regions but we 
were not able to register fluorescence from v' > 26 to the 
X 1 E + state, most likely due to possible predissociation 
or unfavorable transition probabilities. 



B. The a 3 E+ state 



C. Data sets 

Altogether about 6600 transitions to the X 1 E + state 
(130 of them belong to 6 Li 133 Cs) and 180 transitions to 
the a 3 E + state in 7 Li 133 Cs were assigned (see Tables 1 
and 2 of the supplementary materials [3l[). The corre- 
sponding distribution of vibrational and rotational quan- 
tum numbers is shown in Figure [5] For the X^^E" 1 " state 
we observe transitions to v'L = — 50 in 7 Li 133 Cs. The 
set of a 3 E + state vibrational levels covers v'l = 4 — 12 
and the vibrational numbering established in this study 
agrees with that from the theoretical potential. The lim- 
ited number of experimental data, however, might lead 
to a revision of this assignment in the future, if additional 
data on the a 3 E + state are collected. Although the data 
set for the a 3 E + state may seem fragmentary compared 
to our similar studies in other molecules, the collection 
of these data is extremely valuable, first, due to the lim- 
ited possibilities for exciting triplet states in LiCs and, 
second, since a proper description of the Li(2s)+Cs(6s) 
asymptote is only possible if both the X 1 E + state and 
a 3 E + state are treated in a coupled channels manner as 
described in the next section. 



The triplet transitions are easily distinguished by their 
hyperfine structure (HFS) which at our resolution con- 
sists of three lines split by approximately 0.3 cm -1 . Cal- 
culations similar to those in Ref. [35| show that the ob- 
served splitting is well reproduced by the Fermi contact 
interaction model applying the atomic HFS constants for 
7 Li and 133 Cs [36]. In Fig. [2] a transition to the triplet 
state level t>" = 9, A" = 17 [45j is shown together with 
the prediction of the splitting modeled by a coupled chan- 
nels calculation as described in Section HVl A larger por- 
tion of the same progression originating from the v' = 25, 
J' = 17 B X II state level excited by a Q-type (AJ = 0) 
transition is shown in Fig. [3] Some of the transitions 
reach near asymptotic levels in both ground states. The 
vertical dashed lines indicate the prediction of the cou- 



IV. CONSTRUCTION OF POTENTIAL 
ENERGY CURVES 

The self-consistent assignment and fitting procedure 
described above give s rise to accurate pointwise short- 
range potentials [32|]. For the long-range part of both 
potentials we use an extension of the form 



U LR (R) = U a 
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Here Uoo is the energy of the atomic asymptote with 
respect to the minimum of the X 1 E + state potential, C§, 
Cs and C\q are the dispersion coefficients and 
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FIG. 3: (Color online) Portion of recorded progression from the B state level with v' 
transition. Lines not marked on the figure belong to another assigned progression. 
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FIG. 4: (Color online) Observed ro-vibrational levels in the and a 3 E + states of 7 Li 133 Cs and 6 Li 133 Cs. 
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is the frequently applied functional form of the exchange 
energy [37| which is added for the triplet state and sub- 
tracted for the singlet state. 

The pointwise short-range and the long-range poten- 
tials are connected at a point R a ensuring a smooth tran- 
sition between both potential branches. R is chosen as 
described in Refs. [27], [38|, the Cr and Cr coefficients are 
fixed to their theoretical values [H, HO, E[ , 7 and (3 are 
estimated using the ionization potentials for Li and Cs 
[42| according to Ref. [37| . while Uoc, C10 and A ex are 
adjusted during the fitting procedure. 

In the first step of the fitting procedure we adjust only 
the pointwise part of the potential for the X 1 E + state. 
The experimental transition frequencies are fitted by ad- 
justing the parameters of the pointwise potential and the 
term energies of the excited levels. 

As we choose the origin of the potential energy at the 
minimum of the X^" 1 " state potential, in the second step 
of the fitting procedure we need to determine the term 
energies of the a 3 E+ state levels with respect to the X 1 ^" 1 " 
state. We do this using spectra where we observe simulta- 
neously progressions to the singlet and the triplet states 
originating from a common upper state level. For a given 
singlet ground state potential we calculate the term en- 
ergies of the triplet state using the energy of the excited 
level, determined in the previous step, and the progres- 
sion to the a 3 E + state from this level. These a 3 E + term 
energies are then used in order to fit the potential param- 
eters of the triplet state. In this way we ensure always 
a proper position of the a 3 E + state with respect to the 
X^ 4 " state. 

In order to treat the hyperfine structure of the spectral 
lines of the triplet state we checked the experimental data 
that the splitting within our resolution is independent of 
the vibrational and rotational quantum numbers v" and 
N" (except for very few cases, which we discuss below). 
Therefore, we use the central component of the struc- 
ture for identification of the transition (see Fig. [2|) . We 
convert the observed frequency to term energy, and take 
into account the shift (-0.035 cm -1 , which means the hy- 
perfine level is more deeply bound than the unperturbed 
level) of the selected hyperfine component from the un- 
perturbed, hyperfine structure free one, thus we fit with 
the constructed term values a Born-Oppenheimer poten- 
tial. 

For high vibrational levels, and especially in cases of 
close approach of singlet and triplet levels with the same 
rotational quantum numbers, significant deviations from 
the Born-Oppenheimer picture can be expected [27l. l28l|. 
In 7 Li 133 Cs this is most pronounced for v'x = 49 in the 
singlet and u" = 10 in the triplet state (see Fig. [3j . The 
deviations of the transitions determined by the Born- 
Oppenheimer potentials (indicated with thick blue bars 
in Fig. [5]) from the experimental ones are significant and 
reach 0.06 cm" 1 for v" x = 49 and J' = 24. Fig. [5] shows 
the development from low J to high J. 

Therefore, in a third step, the Born-Oppenheimer 
potentials are refined by extending the single chan- 



nel approach from the first two steps and applying a 
coupled channels calculation as discussed in detail in 
Refs. [13, HH . Briefly, we calculate the difference be- 
tween the single channel and coupled channels eigenval- 
ues and subtract these differences from the experimen- 
tally observed transition frequencies in order to obtain 
the frequencies which would be observed without hyper- 
fine coupling between the singlet and the triplet ground 
states. Next, these frequencies are used in a combined 
fit to adjust the parameters of the long-range extensions 
of the potentials as well as in separate fits (described for 
the first two steps) of the short range pointwise part of 
the triplet and singlet potentials. Finally, the whole fit- 
ting procedure is repeated until the frequencies predicted 
with the coupled channels calculations agree with the ex- 
perimental observations. The importance of the coupled 
channels calculation is illustrated in Fig. [5j where the 
predictions for the coupled system are shown with thin 
red bars together with the experimentally observed lines 
and the single channel predictions (thick blue bars). The 
degree of mixing between the triplet and singlet levels is 
characterized by the expectation value of the total spin 
operator \S\ which is indicated for each line in Fig. [5) 

In Tables U and |TT] the fitted potential energy curves of 
the X 1 E + and a 3 E + states are given. The dispersion co- 
efficients Cq and are taken from Refs. [3(| |4l[. With 
the present data sets we are also able to reproduce the 
experimental data with the same quality of the fit by fix- 
ing these coefficients to the values from Ref. [39]. The 
reason for choosing the more recent values is that in this 
case the fitted C10 coefficient differs from the theoretical 
prediction by only -17 %, whereas if the leading disper- 
sion coefficients are fixed to the values from Ref. [39j the 
difference reaches +68 % (the derived C%o amounts in 
this case to 2.39 x 10 10 cm -1 A 10 ). The selected set gives 
good consistency with the expected accuracy of most re- 
cent calculations of dispersion coefficients. 

The potential curve at any point R < R is defined 
by the natural cubic spline function through all points 
listed in Tables Q] and [II] For R > R D the long range 
parameters and expressions (P) and @ should be used. 
For convenience of the reader, we give in Tables [J and 
UT1 also T e (energy of potential minimum), R e (equilib- 
rium distance) , D e (dissociation energy) and the dissoci- 
ation energy Dq with respect to the lowest rovibrational 
level (v" = 0, J" — 0). While the model parameters for 
the potential are listed with all relevant figures neces- 
sary to reproduce the model with sufficient precision, for 
the dissociation energies as physical quantities uncertain- 
ties have been estimated according to the data situation. 
Especially for the a 3 E + state the uncertainty of the dis- 
sociation energy is fairly large because of only few data, 
especially no data are available yet for the lowest vibra- 
tional levels. 

The derived X 1 E + state potential describes the experi- 
mental transition frequencies involving 2400 energy levels 
of the ground state with a standard deviation of 0.0057 
cm" 1 and a dimensionless standard deviation of a = 0.62. 
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FIG. 5: (Color online) Observed transitions to v'x = 49 and u" = 10 for different J' as indicated on the figure and excited 
by a Q-type transition. An extended portion of the spectrum for J' — 17 for several w" is shown in Fig. [3] For increasing 
J' the singlet and triplet lines approach each other, resulting in strong perturbations of the line structure and the intensity 
distribution within the HFS of the triplet line. Thick blue bars indicate the line positions as predicted by the Born-Oppenheimer 
potentials. Thin red bars indicate the prediction of the coupled channels calculation. The numbers close to these bars are the 
corresponding expectation values of the total spin operator \S\. For convenience spectra are adjusted such that the largest 
peaks have similar intensities. 



The high standard deviation of the fit compared to the es- 
timated error limits of 0.003 to 0.005 cm -1 from the typi- 
cally applied experimental resolution of 0.03 - 0.05 cm -1 
arises from a relatively large number of supplementary 
spectra (giving rise to about 26 % of the identified tran- 
sitions) recorded at lower resolution (0.1 cm -1 ) which 
are also included in the data analysis. The standard de- 
viation for the experimental data with uncertainties less 
than 0.005 cm -1 (about 3600 transitions) amounts to 
0.0030 cm -1 and a for this case is 0.88. The increase of 
the dimensionless standard deviation is most likely due 
to overestimated error limits of the low resolution lines. 
The quality of the triplet state potential is assessed by 
comparing the experimental term energies with the cal- 
culated eigenvalues. The standard deviation amounts to 
0.0044 cm -1 and the dimensionless standard deviation is 
0.51. 



In addition to the potential energy curves a set of Dun- 
ham coefficients was fitted to the data for the singlet 
ground state. These coefficients are given in Table 5 of 
the supplementary materials [3l| and describe the exper- 
imental data for all rotational quantum numbers and a 
reduced set of v x : < v' x < 45. 



V. CONCLUSION 

Highly accurate potentials for the singlet and triplet 
ground state were derived, from which one can read off 
the quality of the ab initio result [26[. Because a graph- 
ical comparison is generally too rough we compare in- 
stead two quantities, namely the dissociation energy D e 
and the equilibrium internuclear separation R e . For the 
X 1 ^ state we find D e =5875.455 cm" 1 and i? e =3.6681 A 
with the corresponding ab initio results being 5996 cm -1 
and 3.615 A, respectively [26j]. For the a 3 E + state the 
present work reveals D e =309cnT - 1 and i? e =5.2472A, 
while the corresponding ab initio calculations give, re- 
spectively, 307 cm" 1 and 5.229 A. 

The energy difference in the case of the singlet state 
could correspond to at least one vibrational level more 
than the total number of levels accommodated in the 
derived potential (55 for J=0), but for the triplet state 
the agreement is surprisingly good. Such precision is 
rather good and helpful for guiding the spectroscopic 
assignment. The amplitude of the exchange interaction 
which can be estimated from the difference of the the- 
oretical asymptotic singlet and triplet potentials comes 
close (30%) to the value from the fit. 

Cold collisions were studied for Li+Cs pairs through 
sympathetic cooling by Mudrich et at 21| and trap loss 
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TABLE I: Pointwise representation of the potential energy 
curve for the X E + state of LiCs. See also Table 5 of the 
supplementary materials [31J. 
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measurements by Schloder et al. [201 ]. The latter work 
gives loss rates for processes where excited states are in- 
volved and thus cannot be related directly to cold colli- 
sion calculations which are now possible with the ground 
state potentials reported in this paper. The former work 
derives the cross section for elastic scattering of unpolar- 
ized atomic pairs in the hyperfine ground states fu = 1 
and fcs = 3 to be 8(4) x 10 -12 cm 2 assuming Wigner's 
threshold behavior for s-wave scattering, i.e., the cross 
section is independent of collision energy. In this collision 
process the channels /=2, 3, and 4 of the total atomic 
angular momentum / = fi^ + fc s are involved. Using 



TABLE II: Pointwise representation of the potential energy 
curve for the a 3 E + state of LiCs. See also Table 5 of the 
supplementary materials [3lJ. 
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the potentials reported in this work the cross sections 
of these elastic channels were calculated for an energy 
range up to 100^K. The values for / = 3 and / = 2 
are at least an order of magnitude smaller than that of 
/ = 4, thus only this channel should be taken into ac- 
count for the comparison to the experimentally derived 
value which has 50 % uncertainty. The value of this cross 
section varies only from 2. 5 to 2. Ox 10~ 12 cm 2 from zero 
to 100/iK energy. Thus the threshold law is sufficiently 
well fulfilled. By weighting the cross section of / = 4 
with the statistical weight 9/(9+7+5) according to all 
existing channels we get as cross section of the unpolar- 
ized collision 1.1 x 10 -12 cm 2 , which is about a factor 8 
smaller than the one derived from the experiment, but 
within two times the given error. The potential for the 
singlet ground state is well determined by a large body 
of data (2397 levels), but the triplet ground state was 
only determined by 89 levels. Thus we believe that the 
difference is not directing to a discrepancy between both 
results but ask for more spectroscopic data or more pre- 
cise cross section measurements or a direct observation 
of Feshbach resonances which can be incorporated in the 
fit of potential functions to obtain a full description of 
the spectroscopy and the cold collisions of Li + Cs atom 
pairs. Calculations of scattering length for the singlet 
and triplet states show that for the singlet ground state 
it is already well determined by the present study and 
will be 50(20) a (atomic unit a = 0.529 x 10~ 10 m), 
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but in the case of the triplet state large values with dif- 
ferent signs were obtained during the evaluation with po- 
tentials represented with spline coefficients or with piece- 
wise analytic functions as used in some of our other work 
(see e.g. [HI). Thus, such values are not yet reliable 
and calculations of Feshbach resonances would be of no 
value. To make such calculations reliable more spectro- 
scopic data for the low vibrational levels of the a 3 E + 
state and of near asymptotic levels of both ground states 
is needed. In addition to collecting data by exciting to 
high lying levels of the B 1 ]! state we searched for, but did 
not yet find, other excitation channels which will lead to 
combined fluorescence to both ground states, especially 
to asymptotic levels. For improving predictions of ap- 
propriate excitations we started new experiments to get 



precise data of various excited states. The analysis of 
the states B 1 !! and D 1 !! is almost complete and will be 
published in a forthcoming paper (34[. 
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